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Abstract 

In this paper we statistically analyze the Fokker-Planck (FP) equation of Schramm-Loewner 
evolution (SLE) and its variant SLE(k, p c ). After exploring the derivation and the properties of 
the Langevin equation of the tip of the SLE trace, we obtain the long and short time behaviors 
of the chordal SLE traces. We analyze the solutions of the FP and the corresponding Langevin 
equations and connect it to the conformal field theory (CFT) and present some exact results. We 
find the perturbative FP equation of the SLE(re, p c ) traces and show that it is related to the higher 
order correlation functions. Using the Langevin equation we find the long-time behaviors in this 
case. The CFT correspondence of this case is established and some exact results are presented. 


1 Introduction 

Schramm-Loewner evolution (SLE) is a stochastic process which aims to classify the two-dimensional 
critical models by focusing on their geometrical quantities mm- The idea of classification of probability 
measures of random curves growing in the upper half complex plane was first proposed by O. Schramm 
Q] . This technique has provided us with a new interpretation of the traditional conformal field theory 
(CFT) and Coulomb gas approaches. In Conformal field theories (CFT) the emphasis is on the algebraic 
properties of the critical models dictated by their conformal invariance [3]. Therefore it seems reasonable 
that there is a connection between this theory and SLE. This relation was found in the pioneering work 
of M. Bauer and D. Bernard [1] in which it was found a relation between the central charge of CFT 
and the diffusivity parameter («) of SLE. 

Many papers have been published concerning the analytical and numerical properties of interfaces of 
various two-dimensional (2D) critical models described by SLE theory |l* 5] (i [7]. The statistical 
analysis of the various aspects of these critical curves contains valuable information of the model in 
hand. In the numerical studies of the less-known 2D critical models, the main concern is to extract the 
diffusivity parameter (referred to as k) of their critical interfaces to relate the underlying model to a 
CFT class 0 0. Therefore the dependence of statistical observables to k is of paramount importance 
in this area EMU- As an example of these statistical observables, one can mention Fokker-Planck (FP) 
equation (18] , left passage probability }9], crossing probability [2], etc. which are directly related to the 
diffusivity parameter (k) of SLE. In spite of many analysis addressing the FP equation for the SLE 
traces in various papers mm, a complete and comprehensive analysis of this subject is missing in 
the literature yet, e.g. the actual behavior of the solutions of this equation in the whole domain of 
the model and a comprehensive investigation of the FP equation for the dipolar and radial SLE and 
SLE(k, p c ) E The question what the statistical properties of the tip of the growing random curves are 
in the SLE(k) and SLE(n, p c ) is the main concern of this paper. The importance of the statistical 
analysis (such as FP equation) of SLE(ac, p c ) can be understood, e.g. for the case in which the curve in 
its growth process experiences some boundary condition changes (bcc) in some point on the boundary 
of the domain of growth. This situation which occurs in most cases in the numerical analysis of the 

1 It is customary to use the phrase SLE(re,p). In this paper we use SLE(k ,p c ) in order not to be confused with the 
probability density p{x. y, t) to be defined in SEC [ 3 ] Our analysis is limited to the case p c = k — 6. 
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less-known models should be investigated within SLE(k, p c ) theory. Such analysis may help to analyze 
these behaviors. 

In this paper we analyze the various aspects of the FP equation of SLE traces. To this end, we 
firstly show how to employ the backward SLE equation and find the Langevin equation for the spatial 
coordinates of the tip of the SLE trace. We analyze these equations in short and long time limits and 
show that for ordinary (chordal) SLE in long time both horizontal and vertical components of the tip 
of the SLE trace behave as 1 2 (t = time) as expected. To do so, we find the behaviors of SLE traces for 
long times in two ways; directly solving the FP equation and using the Langevin equation. We then 
connect the FP equation to the CFT and find some exact behaviors. For SLE(/c, p c ), however there 
is other preferred point on the real axis, causing the situation be different. To analyze this case we 
first find the behavior of the curves for short times perturbatively. The long-time behaviors are also 
investigated by analyzing the Langevin equation. We finally report some exact behaviors of this case 
by means of its CFT counterpart. 

The paper has been organized as follows: in SEC [2] we present a short introduction to SLE(k) and its 
variant SLE(k, p c ). The SEC [3] has been devoted to the derivation of FP equations for the SLE traces. 
FP equation for the SLE(k,/3 c ) is derived in SEC [4] and the short and long time limits of this equation 
is derived in this section. 


2 SLE 

As a well-known fact, the critical 2D models have special algebraic and geometrical properties. The 
algebraic properties of these models are described within the conformal held theories. However these 
theories are unable to uncover the geometrical features of these models since it concerns the local 
fields defined in these models. According to Schramm’s idea one can describe the interfaces of two- 
dimensional (2D) critical models via growth processes named as SLE. Thanks to this theory however, 
a deep connection between their local properties and their global (geometrical) features has been dis¬ 
covered. This theory aims to describe the 2D critical models by focusing their interfaces. These 
non-intersecting domain-walls are assumed to have two essential properties, conformal invariance and 
the domain Markov property. In the SLE theory one replaces the critical curve by a dynamical one. 
We consider the model on the upper half plane, i.e. H = {2gC,Sz>0}. Let us denote the curve 
up to time t as 7 1 and the hull K t as the set of points which are located exactly on the 7 1 trace, 
or are disconnected from the infinity by 7 1 . The complement of K t is H t := H\K t which is simply- 
connected. According to Riemann mapping m theorem there is always a conformal mapping gt{z) 
(in two dimensions) which maps H t —> H. The map gt(z) (commonly named as uniformizing map, 
meaning that it does uniformize the "u trace to the real axis) is the unique conformal mapping with 
gt(z) = z + — + 0{\) as 2 —>■ 00 known as hydrodynamical normalization. Loewner showed that this 
mapping satisfies the following equation 

UEHH]: 

a,9,(2) = (1) 

with the initial condition gt(z) = z and for which the tip of the curve (up to time t) is mapped to the 
point on the real axis. For fixed z, g t {z ) is well-defined up to time t z for which gt{z) = Ct- The more 
formal definition of hull is therefore I\ t = {z £ H : t z < t}. For more information see [2, 0; :IU H5] . 
For the critical models, it has been shown pj that (referred to as the driving function) is a real¬ 
valued function proportional to the one-dimensional Brownian motion = y/TiBt in which n is known 
as the diffusivity parameter. SLE aims to analyze these critical curves by classifying them to the one- 
parameter classes represented by k. 

There are three phases for these curves: 

Dilute phase: for 2 < k < 4 the trace is non-self-intersecting and it does not hit the real axis. 
Therefore in this case the hull and the trace are identical: I\ t = 7 t . 
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Figure 1: Schematic plot of the SLE trace in presence of a bcc point x°° (boundary change has been 
shown Schematically by (+) and (—)). 



Figure 2: Schematic plot of the SLE trace which ends on the real axis at x°°. 


Dense phase: for 4 < k < 8, the trace touches itself and the real axis so that a typical point is surely 
swallowed as t —» oo and I\t ^ "ft- 

Space-filling phase: for k > 8 the curve is space-filling. 

There is a connection between the two phases: in the dense phase the frontier of K t (the boundary of 
H t minus any portions of the real axis) is a SLE curve. It has been shown j5] that this simple curve is 
locally described by the diffusivity parameter k = ^ in the dilute phase. 

The correspondence to CFT is via a simple relation between diffusivity parameter n and the CFT- 
central charge c, i.e. c = ( 3k - s )( 6 - k '> g]. 


2.1 SLE(ft, p c ) 


Let us now suppose that there is a preferred point on the real axis, namely at x°° as indicated in FIG 
[l]and[2j We are going to analyze the random curve in presence of this point. Note that the conformal 
invariance of the curve has been explicitly broken, however the curve is yet self-similar. SLE(k, p c ) is 
a variant of ordinary (chordal) SLE which analyzes such curves. We consider the upper half plane as 
above. In this theory the parameter k is the same as above. The parameter p c has to do with the 
boundary conditions (be) imposed. The point Xoo on the real axis mark the location where boundary 
condition changes (we name it as Xoo since this point may be considered as the point at which the 
random curve returns as t —> oo, see below). The exact value of p c depends on the boundary change at 
this preferred point. In the corresponding CFT, the bcc imposed in Xoo corresponds to putting a bcc 
operator with conformal weight h Pc = at that point f\mm- 

For these curves the stochastic equation is the same as formula (!]). The driving function acquires 
however a drift term with respect to the chordal SLE: 


d( t — \fndBt + 


Ct ~ x' 


-dt. 


( 2 ) 


2 Another bcc operator should be put at the origin with the conformal weight h 0 = K 
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Figure 3: Schematic representation of the conformal map ht(z) and the backward conformal map 
ft{w) — Ct on the SLE trace. 


in which B t is one-dimensional Brownian motion and x%° = <?t(:Coo) is the transformation of the point 
a,’oo under the action of g t . This situation can be generalized to a set up with the arbitrary number of 
bee’s. 

Real-axis to real-axis random curves: 

As an example of the above discussion, we consider the FIG [2] in which the mentioned preferred point 
is the ending point of the curve (which does not go to the infinity, but after the exploration process 
turns to the real axis). For this set-up, by using the conformal map <fi = x ao z/{x^ — z), one can send 
the end point of the curve to the infinity. Since our model is conformal invariant, we can deduce that 
the function ht = 0og t o<(> _1 describes ordinary (chordal) SLE from which one easily concludes that gt 
satisfies d t g t = 2/({<p\g t )((f>(gt) — Ct)}) in which (j>'{z) = The other criteria is to make the mapping 
hydrodynamically normalized. It has been shown HU that the stochastic equation of hydrodunamically 
normalized gt is the same as Eq. (|T]) with the price that Q acquires a drift term as follows HZ]: 

d( t = \fndBt + —^^— rdt (3) 

Ct - gt{xoo) 

In the other words, this stochastic function is the driving function of the SLE(tc,p c ) with p c = k — 6. 
Thus for the critical curves from boundary to boundary, the corresponding driving function acquires a 
drift term. This generalization of SLE can be generalized further to have multiple preferred real axis 
points and multi-critical curves [§J. For review see also references (atm. 

The other example is the traces which terminate surely to the real axis in the interval [xo,oo). For 
these curves of course the point xn is regarded as the preferred point for which p c = [k — 6)/2 yields 

dipolar SLE(/s) 01 El HU- 

3 Derivation of Fokker-Planck Equation 

In this section we are going to find the Langevin equation for the tip of the chordal SLE trace. Let us 
first define the shifted conformal map as follows: 


ht(z) = g t {z) - Ct 


(4) 
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for which one can easily verify that 


h t 1 {w) = 9 t 1 (w + Ct)- 


(5) 


In the above equation = means the equality of the distributions of stochastic processes. The differential 
equation governing ht is: 


d t h t (z) = 


h t (w) 


- dtCt 


( 6 ) 


in which ho(z) = z. One can retrieve the SLE trace by the relation 7 (t) = lim^o {Ct + * e ) = 
linu-^ gf 1 {z) and find the differential equation of the tip of SLE trace. The equation governing g t _1 
is very difficult to solve. There is a way out of this problem, using the backward SLE equation. It is 

defined as follows (note that Ct = C-t = —C t) HHj : 


dtft{w) = -- 


ft (w) - C t 

It has been shown that the probability distribution of ft is the same as g~ t 1 [HI, i.e. 


ft(w) - c t = 9t 1 ( w + C t) = h t 1 (w). 


(7) 


( 8 ) 


The schematic representation of this equation has been sketched in FIG [3] Therefore the tip of the SLE 
trance can be obtained 7 t = lim^^o gf 1 ^’ + Ct) = lim^^o hf 1 (w). Now we can find the trajectory of 
the tip of the SLE trace using the backward equation 0 for z t = x t + yt which is (notice that Re( 7 t ) 
and Im( 7 t) have the same joint distribution as x t and yt) 


d t x t = - - dtCt 


IIt 


a„. 2 y t 

UtVt 2 1 2 

** + Vt 


(9) 


conditioned to have the initial values Xq = u and yo = v in which w = u + iv is the initial point of 
the flow. The Fokker-Planck equation for SLE governs the probability density p(cc, y, t) for the critical 
curves which is the probability that the tip of the curve be at the point (x, y) on the upper half plane 
at the time t. It is defined as p(x,y,t) = (S(x — x t )5(y — yt))- To obtain the FP equation, we notice 
that p{x + dx, y + dy , t + dt) = p(x, y , t). Using the Langevin equations Q of the growing SLE curve 
and after performing some Ito calculations, noticing that for the chordal SLE, the driving function is 
proportional to the Brownian motion (which implies that ((dft) 2 ) = ndt), it can be easily proved that 

m- 

2 (x - Co) Q 2 y 


d t p(x,y,t) = 


:di 


{x - Co) 2 + y 


(z- Co) 2 


y 


2 d V 


p{x,y,t) 


( 10 ) 


in which Co is the point at which the curve starts to grow. This equation has the symmetry (x — Co) ► 
X(x — Co) , y —■► Ay and t —> A 2 t in which A is some (non-zero) positive real parameter. One can see 
this symmetry in the other way. Let us observe the effect of the transformation x —> Xx , y —> Xy 
(setting Co = 0 ) on p(x,y,t) which yields (S(Xx - x t )6(Xy - y t )) = -h{5{x - \x X 2 T )5(y - xy^r)) in 


which T = A 2 t. Noting that t£a 2 t = xt and yj/A 2 r = Vt (see Eq. 

p{Xx,Xy,X 2 t ) = j^p(x,y,t). 


9), we find that 


Considering also the translational symmetry of the problem, i.e. x —> x + a and Co 
its general form to be: 

p(x, y, t) = r~ 2 g (xt, v) = t~ x f (\t, v) 


( 11 ) 

Co + a, we expect 
( 12 ) 
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in which \t — 
determined by 
[lO] we see that 


(x-Cof+y 


ry = a, r = yj{x — Co) 2 + V 2 and / and g are some functions to be 

Eq [lO] (not to be confused with mapping f t defined before). By setting Eq 12 into Eq 
it is reduced to the following form: 


h xx d xf + 2h xv d x d vf + h m d vf + h x d xf + h v d vf + / = 0 


^XX 


2k 


V 2 X 

1 + ry 2 


h 

h 


_ k (1 + ry 2 ) 

2 x 

xv = 


h x 


K + X + 


4 (ry 2 — 1) 
r? 2 + 1 


hrq 



(13) 


It can easily be checked that this equation is a hyperbolic PDE, i.e. 
£ = ^ (so that x - £ 0 = vJ § and y = J±), and f(x,v) = 


h £ — h xx h vr) = 0. By defining 
F(£, ry) we obtain the canonical 


form: 


1 

2 


K0ri F 



1 + rj 2 . ?y 2 - 1 

-r-b 4—^-— 

£ ry 2 + 1 



5 ? E + 


4?y£ 

1 + ry 2 


d v F + F = 0 


(14) 


Before investigating this equation, let us explore the properties of the regime of our interest, i.e. long 
(non-primitive) times. In this limit we expect that Xt yt- Consequently we are interested in the 
large ry limit of the above equation. In this limit we have: 


Ud^F - d^F +-d v F + y = 0 (15) 

z ry £ 

To solve this equation we consider the factored form of solution, i.e. F = II(ry)E(£) from which we 
conclude that: 


f - (A + |)E = 0 

\ |/c5gn+^n-An = o 

In this equation A is an arbitrary real number. The solution of the above equations are 


(16) 


E(£) = So ( | 


oKZSo) 


Il(ry) = ry 2k I k-s [ i\j —— ry ] + ry 2 » K K -s 



(17) 


in which I and K are modified Bessel functions of first and second kind respectively. The expected 
behavior of these functions (being real in each interval) is satisfied only when A = CP which yields 
n(ry) = ry s - and causes £(£) to have a maximum at some point and uniformly go to zero at large 
values of £. Therefore for long times we obtain that 


P(x,y,t) 


po m) 

t 77 ( 8 — K )/( K ') 


Po 

t 


8 -«)/(«) 


m) 


(18) 


in which p 0 is the proportionality constant. This result reveals some important behaviors. For k < 8 
the probability of touching the real axis boundary, i.e. ry = | — 1 ► ±00 goes to zero as , but for 
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k = 3 



Figure 4: The solution of the equation pFTT| for the SLE curve with diffusivity parameter k = 3. In these 
graphs F has been plotted in terms of % (77 = 0). 


k = 8 there is a finite (maximum) probability of touching the boundary which is an important effect 
in SLE, having its root in the fact that for this case the curve is space-filling. For large amounts of 77 
also, p(x,y,t) ~ lim A _> 0 - £e A ^ goes to zero for large £ values. 

To be self-contained, we have numerically solved the equation [Ll] To this end we have used the finite 
difference method. We have considered a lattice with dynamic mesh sizes (bounded to 0 < y < 10 3 
and —2 x 10 3 < 77 < 2 x 10 3 ). Due to the singular behavior of Eq. (13 1 near y = 0 and 77 = 0, one 
should mesh the space near these lines more exact (smaller) than the other points [^j The results have 
been indicated in Figs. BE As is seen in Fig. [4j when 77 = 0 (x — > 0 and t = finite), F has a single 
peak at some y. In fact we have observed that these peaks occur at yo = x °+ Vo ~ ^P- = a(k) in which 
a(/v) ~ (k) - 5 . This fact shows that for larger amounts of k. at a fixed time, the vertical growth (along 
the y axis) of the growing curve decreases. Figure [ 5 ] shows this dependence for large 77, i.e. |x| y 
and t = finite. The figure shows a very narrow peak appears around Xo — = 6(k) in which &(/?.) ~ k 

which is expected, since for long times Xt — (Eq. [ 9 ]). The symmetry of F with respect to 77 is 

evident in Figs. [6] and [7] in which the behavior of F has been sketched in terms of 77 for the cases y = 1 
and very large amount of y respectively. For y = 1 (r 2 = t) it is shown that F has a peak around 
77 = 0 and decreases (for positive 77 axis) as 77 increases and for 77 —> 00 (2000 in the graph) it tends to 
zero according to the relation obtained above. For y > 1 (1000 in the Fig. [7]) the sensitivity of F to 77 
(= cot 6) becomes negligible and its absolute value becomes very small. 

Before closing this sub-section, we notice that from the general solution of F(£,? 7 ), we can also 
calculate (x n ) and (y n ) as follows: 


/ OO /»00 

dx / dyx n p(x, y, t) 

-00 J 0 

=¥ 11 - 1 ?^ 

/ oo /»oo 

dx / dyy n p(x, y, t) 

-00 J 0 


1 n 

= -£2 

2 






3 The MATLAB PDE toolbox has been used to solve this equation. 


(19) 


( 20 ) 
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r) = 2000 


8 . 10 . 
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Figure 5: The solution of the equation [13| for the SLE curve with diffusivity parameter k = 3. In these 
graphs F has been plotted in terms of \ (v = 1000). 



r) t3 

Figure 6 : The solution of the equation [T3| for the SLE curve with diffusivity parameter k = 3. In these 
graphs F has been plotted in terms of 77 (% = 1). 



Figure 7: The solution of the equation |13| for the SLE curve with diffusivity parameter k = 3. In these 
graphs F has been plotted in terms of 77 (% = 100). 




































































(x n y m ) (t) 


d.T / dyx n y m p(x,y,t) 

> Jo 

V n F(£,v) 


-i /*oo /»oo 

f i(n+ m ) / dr? / d £ 

" J—oo J 0 


0 ^i(n+m )+2 


(21) 


We see that the n-th moment of displacement scales as ( r n ) ~ t 2 . We can learn more about F from 
these relations. Since Eq[l4]is invariant under y —► —y, one can conclude that its solutions are classified 
into even and odd solutions with respect to the ?y component. By setting n = 1 in Eq |19[ and using the 
fact that (x) = 0 we can easily find that only even solutions should be taken into account, i.e. F(£,y) 
is even in y component (note that we can not directly generalize this conclusion to the limiting case 
F ~ i 7 - ( 8-k )/( k ). In fact as we see later this function should be of the form (1 + ?y 2 ) _ ( 8_K )/( 2K ) which 
is even function of y and for large amounts of y becomes the above result). 

These results are valid within the Fokker-Planck equation and should be verified by the Langevin 
equations [9] which is the subject of the next sub-section. 


3.1 Langeving equation for the long-time limit 

In this sub-section we obtain the behavior of x(t) = \J(x 2 ) (t) and (y) ( t ) for long times, using the 
Langevin equation ( Eq. |9j). In this limit we have x *^_ y 2 —> 0 and therefore x t ~ Therefore 

x(t) — in accordance with Eq. 19 To find the behavior of yt, we notice that in long-time limit, 


the probability distribution for x t = y/KB t is P{x,t) = 
replacing the r.h.s with its average, we obtain: 

1 r 


= e Using the second line of 

I 0 


Eqg 


and 


d t y t ^ 


2 Vt 


re v^Ax 


V2t TKt 7-oo Vt + x 2 
=>e _x * ({Kt)idtiVtUxt)^ ~^J e~ v ~dy 

in which \t = O ne can easily obtain the differential equation governing Xt- i.e. 

' d — K 


( 22 ) 


dt (tdtXt) - 2 txt{dtxt) - 


2k 


+ Xt dtxt = 0 


(23) 


If we represent this function again as Xt = At 2 f(t ) (in which f(t) is proportional to yt), then the 
following equation is obtained (/' = ^/(i) and f" = ^sf(t)): 


It can easily be checked that the answer of this equation is /(f) ~ t 2 . The conclusion is that in the 
long-time limit we have: 

f x t = yfnB t ~ fs 
\ yt ~ t 2 

which is in agreement with the results of the previous sub-section. 


(25) 


3.2 CFT background 

We have found the solution of the FPE of SLE growing curves in long-time limit in the previous sections. 
There is another way for exploring this problem within conformal filed theory (CFT) techniques which 
is the subject of this section. 
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Let us define the function P(x, y) as the probability distribution of curves passing the region (x, x + dx) 
and (y, y + dy). It is defined as: 

poo pOO 

P{x,y) = p(x,y,t)dt = (6(x - x t )S(y - y t )) dt (26) 

Jo Jo 

Throughout of this sub-section we name the point at which critical curve starts to grow as C 0 which is 
important in our analysis. The behavior of this function under the transformation (x — (o) —t A(x — Co) 
, z/ —J" Aj/ is obvious: 


pOO 

P(x>y)\(x-0->\(x-Q,y->\v —*■ / p(Ax,Xy,t)dt 

JO 

1 f°° 

= J V ’ T ) x a:dT = p ( x > y) 

and its more explicit form is: 

pOO 

P{x,y)= / t~ 1 F{^ t ,y)dt 
Jo 


(27) 


F&V) 


(28) 


d C 


showing that it has zero dimension. 

Let us investigate the properties of P(x, y, e) which is the probability that the distance of the point 
z = x + iy from the SLE curve be e. In the other words, it is the probability that the curve passes 
a disc with radius e centered at z = x + iy. This probability is compatible with P(x, y) with a finite 
area (disc with radius e) around. From the CFT point of view and noting that the original model 
has the conformal symmetry, one can easily understand that P(x,y,e) is proportional to the following 
multi-point correlation function: 


P(z,z,(o,x°°) 


(f(x, y)F{x , -yMCoMz 00 )) 

(tHCoMz 00 )) 


(29) 


In this equation, F{x, y) is a detector operator whose value is non-zero in a field configuration only if 
the mentioned curve (defined in the background of other fields) passes the infinitesimal region around 
the point (x, y) and x°° is the point on the boundary at which the SLE curve ends (in other words, we 
have two boundary changes, one at Co and the other at x°°). The operator F(x, —y) is the image of 
F[x, y) with respect to the real axis whose presence in the correlators in the boundary CFTs (BCFTs) 
is customary, ip is the boundary changing operator for the CFT corresponding to underling SLE whose 
conformal weight is /ii(ft) = with second level null vector [B]: 


(^L\-2L_ 2 )iP = d 


(30) 


in which L n ’s are the generators of the Virasoro algebra satisfying 


[J^rif I'm] — 'Wl)J J n+rn “1“ l)^n+m,0* 


(31) 


and c is the central charge of the underlying CFT. Using Ward identities, Eq. 30 leads to the following 
equation for the mentioned correlation function [3j: 


(^C\ - 2 £_ 2 ) (F(z, y)F(x, -yMCoM* 00 )) = 0 ( 32) 
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with 


C -1 = 


*-» = £ 


d(o 


hi 


d 


{zi- Co) 2 Zi-CodziJ' 


(33) 


in which the sum is over each field in the correlation function Eq. (29) except Co- From global conformal 
invariance we know that 


Q(z,z, Co, x°°) = (F(x,y)F(x,-y)^(Co)ip{x°°)) = 


1 


so that Q satisfies the equation 


*d 2 Co - 2h F Re | - 


1 d 


1 d 


(Co - z°°) 


d 


2/h(k) ■^ > ( 2 ’ z >Go,x ) 


(34) 


/ii(re) 


\(z-Co)*J z-( 0 dz z-(odz x°° - Co dx°° (x°° - Co) 2 . 


Q = 0. (35) 


In this equation h F is conformal weight of F. From general grounds of CFT, we can easily find that 
the four-point correlation function in Eq[29]has the following form: 

p = y-Jfcjr + f fclW^OO _ ^ 0 )f/»i(«)+S^[(a; - x°°) 2 + y 2 ]^ h F*)-h h F[( x _ ^ 0 )2 + y 2 ]^h 

1 


y 


2hi 


M u ) 


(■g—Co)(-g—a 


(36) 

(?/ = lim^oo^oo it) and /io and /i are some functions to be determined. 

(37) 


in which u = Re-=- r 

y(x“-Co) 

For the limit x°° —> oo the above equation reduces to the following form: 


-dl 0 -2h F l\e( , 1 ^' 0 )+—L-lL' + -l^£- Q = 0 . 


V(^-Co) 2 / z-Codz z-( 0 dz_ 


and Re[u] = rj and P(z , z, Co) ~ Q(z, z, Co) = 2 k F h(rf). In this case, by substituting Eq 36 in Eq 
obtain (Co = 0): 
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we 


K 

2 


2x 

x 2 + y 2 


2 y 

x 2 + y 2 


d y + 


8 h F y 2 
(x 2 + y 2 ) 2 


h{u) = 0 


^ (l + V 2 ) 2 d 2 +4?? (l +r/ 2 ) d v + 8 h F h(rj) = 0 


(38) 


Now we note that the conformal weight of P(x, y, e) should be zero. The only scale which may com- 

/ \ 2 h f 

pensate the factor y~ 2hp in front of P is e so that it should be of the form P(x, y, e) ~ ( ^ j h{rj) = 
ynrfh(r],e) so that 2h F h(rj,e) = ed e h(r],e) which yields: 


^ (1 + rj 2 ) 2 d 2 + 4ry (1 + rj 2 ) d v + 4ed e 


h(u, e) = 0 


(39) 


This equation can be obtained directly by employing the SLE equations, noting that under each confor¬ 
mal transformation gst{z), the mentioned disc is transformed as e —> \g' St {z)\ e ~ (1 — 2<5tRe(l/.z 2 ))e |2J. 
If we use the trial solution h{rj) = (1 + r? 2 ) 7 for the above differential equation, we obtain h F = \ — 
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and 7 = . We know that the fractal dimension of the SLE curve is df 

2hp = 2 — df. Therefore we have: 


P(x,y,e) 


(e\ 2 ~ df 1 

\y) (i + 7? 2 )^ 


e 2 -d f 

y>(8 —K) I K 


yW~ 8) 2 


1 + k/8 [ 21 , yielding 


(40) 


First observe that P ~ r §~~ (sin #) - *- 8- ”^ 8 (l + cot 2 9 ) ^ ^)/ 2h ^ = co t 0 ) , showing that 

as k (and also the fractal dimension) increases, P becomes more long-range. The obvious example 
is k = 0 whose curves are straight lines and its probability distribution is expected to have the form 
P ~ r _1 . The opposite limit is re = 8 whose curves are space-filling which leads to the fact that the 
probability distribution is independent of r. The other important result from this formula is that in the 
limit y —» 0 and x = finite, P = 0 except for the case k = 8. Note that this solution is compatible with 
the one obtained in SEC ^ in which it was shown that for rj 1, 11 ( 77 ) ~ ? 7 - ( 8-K )/ K which is exactly 
compatible with the dependence of P(x,y) on 77 obtained above, i.e. (1 + 77 2 )~( 8_k )/ 2k ~ ^-( 8 -k)/k_ 
This is due to the fact that for large limit of 77 , F is separable into II and E and the integral in Eq [26] 
II can be brought out of the integral, maintaining its dependence. 


4 Fokker-Planck Equation of SLE(«, p c ) 

Now we consider the SLE curves which experience a boundary condition change in some point on the 
real axis ( x°°). Specifically we analyze the curves which goe from origin to a point on the real axis. As 


mentioned in section 2.1 these curves are describe by the framework of SLE(k, k — 6 ) in which we use 
the relation of the driving function Eq [3] . The Cartesian components of the tip of SLE growing curve 
satisfy the following Langevin relations: 


dx t = 9 * 9 dt — dQt 


-2 It 


dy t = 


2yt 


■Vt 


dt 


(41) 


d( t = yJlidB t 


-dt 


C t-x?- 

Using this equation and some calculations as above, one can calculate the change of the function 
p(x,y,t). After some Ito calculations it can be proved that: 


dp 

dt 


-d 2 


2x 


, d x - 


2 y 


2/ 


y 


■2 


P + Pcd x { 


1 


(St - X?) 


S(x - x t )S(y — y t )) 


(42) 


in which the second term in the r.h.s. is due to the existence of the preferred point on the real axis. 


4.1 Perturbative FPE of SLE(k, p c ) 


To Evaluate the Eq. 42 we first consider the case x%° Q i- e. for short times of the evolution. In 
this limit x^° ss x°° and Q <C x°° so that we can treat as the perturbation parameter. Expanding 
the second term in Eq[42]in powers of this parameter we reach the formula: 

1 = 1 oo/ O - x t )5(y - y t )) 

(St - X?>) 


= E 

n—0 
00 

-E 

n—0 


1 


St 

-yOO 


(a; 00 ) 


n+1 


5(x - x t )5(y - y t ) 
((St) n S(x - x t )S(y - y t )) 


( 43 ) 
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To the first order of -%■ we have: 

1 - ((4 - r4* W* - **)<% - »*))■ (44) 

^ {xt°) 

Taking x^° ~ a: 00 and using the equation cc* + £t = — f Q ~^yi ds, we can rewrite the above equation 
as: 

—1 1 f 1 —2x 

1 - <(** + y o ^ y*)) 

= + (^ 00)2 ( X P° + 2 -0 ( 45 ) 

j=</ r 2 4,2 ds 6 (x-x t ) 6 (y-yt)). 

Jo x s + 2/s 


To proceed, we calculate I: 

rt 


1 = 


ds 


lo 


—2 --2<50 r/ - £«)<%' - 2 /s)<J(x - Xi)<5(y - y t )) dx'dy' 

x' + y' 


(46) 


= ds 


x' 2 + y 


~^p 2 (x,y,t;x',y ,s) dx'dy 1 . 


In the above formula, p 2 (x,y,t]x' ,y' ,s) = ( S(x' — x s )5(y' — y s )S(x — Xt)S(y — yt)) is the probability 
density of the tip of the curve being in (x,y) at time t and in ( x',y' ) at time s (s <t). The result is: 


dp 

dt 


+ 




( x °°) 2 


2x 


x 2 + y 1 


- — I d, r - 


_^_ q 

x 2 + y 2 v 


p{x,y,t) 


d x {xp(x,y,t) + 2 ds — -k d x p 2 (x, y, t; x', y', s) dx f 

Jo J x' +y' 


dy' 


(47) 


+ 0{{x°°)- s ). 


The above formula can be generalized to any order of the perturbation parameter . The so¬ 
lutions of the above equation needs the information about p 2 and higher order probability func¬ 
tions which are less important than p for early times. For larger times one also may use the hy¬ 
pothesis p 2 (x,y,t;x',y',s) ~ p(x, y, t)p[x', y', s). Let us write every thing up to O ^(x°°) ^ , i.e. 
p(x,y,x°°,t ) = p(x,y,t) + epi(x,y,t) in which e = Therefore, putting this function into Eq. 
up to first order of e we find that: 


47 


dpi 

dt 


k 2 (,k + 2)x a 
2 x+ x 2 + y 2 x 


2 dy + 4 (( K _ 4 )y 2 + 4a: 2 ) 
x 2 + y z r 4 


Pi(x,y,t) 


(48) 


In this case we also have the global scale transformation symmetry and translational invariance, result¬ 
ing to the same form for pi as SEC|3j i.e. pi = t~ 1 il>{Xi 7 l) f° r V’ being a function to be determined. 
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Straightforward calculations yield: 


ci xx d 2 ip + 2a xv d x d ri ip + a vri d 2 ip + a x d x ip + a^d^ip + aip = 0 


_ o . 11 2 X 

&VV - ^ft Q 

1 + if 


ft 

a m = 2 


a xr] = ktj 


(l+? 7 2 ) 


X 


= « + X + 


2((w + 2)^ - 2) 

rf + 1 


( 49 ) 


a„ = (k + 4) — 
X 

CL = 1 “h 


2x(l + r? 2 ) 

which is parabolic PDE as the Eq|l3j 
as the previous sections): 


((k + 2)r] 2 — 2) 


Since CL-^y-j CL xx CL'qr) 


= 0 and its canonical form is (£ = 


1 


-n£d 2 ip — k + 


1 + r] 2 (K + 2)if-2\{ £ 


2 ' '' V ' £ 

+ (k + 4) ^ d v ip + 
1 + rr 


r/ 2 + 1 

£ 


1 + 77" 


((k + 4) + 477 2 


dt-ip 
ip = 0 


(50) 


2(1 + 77 2 ) 2 

If we repeat the same steps as SEC [ 3 ] we will find that 11(77) ~ t]~^ 8+K ^ K , in which we have separated 
ip = III (t?)S i (£) for large amounts of 77. These results are gathered in the following equation: 


p{x,y,x°°,t) 


yS(0 (11(77)+ elli(77)) 

( a ° + T ai 
£ j.00 ^( 8 +«)/ac 


(51) 


in which ao and ai are some proportionality constants which can not be determined in this approach. 
The other way to obtain the behavior of Xt and yt in terms of t more exactly, is by using the Langevin 
equation Eq. |41| Therefore in the next sub-section we analyze other properties of the growing SLE 
traces in dipolar (boundary to boundary) set up, first by its Langevin equations and secondly by means 
of CFT techniques. 


4.2 Langeving equation for the long-time limit 

The long-time limit of Eq 41 ( 2x f 2 ~ 0 and £ t xf) isFj 

_ x t i Vt 

dxt — dQ — —\ffdB t — ^fdt 
Ct 


Vt - [ ^ 2 PBe BS (x,t)dx 

J- 00 x + Vt 


(52) 


in which Pp, etis {x. t) is the probability distribution of the Bessel process, i.e. the first expression of Eq 
[52} In fact the Fokker-Planck equation for this Langevin equation is readily: 

dtPBess = —p c d x (53) 

4 Throughout this subsection be use p c and assume that p c = k, — 6, although it is true for general p c . 
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Figure 8 : The behavior of the hypergeometric and Laguerre functions for k = 4 and p c = 2. 


2 

By changing variable ( = we obtain 

^ 2ti^~ + ^ dC^^ + = °' (M) 

By the transformation Pge ss = C^ e ~^ w (0 we see that w(Q satisfies the hummer’s equation: 



/ 3 K - 2 p c 
^ 2k 



(55) 


whose solution is: 


w{ C) = ^4iFi 


Pc 3 
k ’ 2 


Pc 

5 

K 


C 


+ 5L 



[C] 


(56) 


in which the first term is the hypergeometric function and the second term is the Laguerre polynomial. 
The second independent solution could have been assumed to be ( 2 ^ 1 F 1 [^, K ~^f c , C] > but the 

Laguerre function is more suitable as becomes clear in the following. In FIG [ 8 ] we have presented the 
plot of the two functions stated above. We see that as ^7 = ( =► 00 , the hypergeometric function goes 
to zero, in contrast to the Laguerre function which remains finite. So we conclude that 5 = 0. Now let 
us calculate ( x ) and J^ y i ^ which are: 


c) = \/2nt 


2 y t 

2 +Vt 


/o°° v / C e ~^ 1 F 1 [i-£,|-£,C] dC 
/“e-CiFr [1 -M-^C] d( 


ft 


0 2 Kt^+y\ 




(57) 


r^iFr [l-*,§-*,C] ^ 
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The solution of the first expression of the above formula (x) is: 


(x) = fy/2nt 

Trisect) r2 2 Fi(|,l — & I — & 1) T(1 + f ) ^(f, f + & | + & 1) 


/ = 


4r(? 


ra-f) 




r(! + f) 


( 58 ) 


which shows that in the long-time limit we have ( x ) ~ fs as the chordal case and the presence of the 
preferred point only changes the proportionality constant of this dependence. To obtain the amount of 
y t , we note that the answer of the second equation in Eq |57| is: 


,°° 1 
lo a C + fr 2 


e- C iF! 


Pc 3 _ Pc 

^ ’ o ’ S 

K, 2 K 


d( = G 


1 ^31 ( 0,5 


23 l 0,0,-i + ^f 


(59) 


in which G™” (bi b^ bq” I 2 ) * s Meijer G-function [20- We have (xt = as above and w = lnt): 


1 dxt 
Xt d w 


2«r(| + f) 


G 23 (o’o-i+^IX?) 


It is difficult to find the general solution of the above differential equation, but if we expand 

——--—-TXT-yr = a o + OL\(xt - 1) + 0((xt ~ l) 2 ) 

Xt ( 2 + 2Kr(i + £y) G 23 (o’o-i + felXt 2 )) 

we will be able to determine the above integral as follows: 


(60) 


(61) 


Xt = 


yt 


ao 


yJlKt 


1 H-I — 1 ( 1 T 2 — ^ In t 


a i 




On 


(62) 


which shows the logarithmic deviation from the ordinary SLE solution, i.e. y^ hordal ~ \j2Kt. ao and 
cti in the above formula are: 


ao = 


2 + 8 


Ot!.\ — 


VI c31 2 

2K r (i + fe) LT23 (,0,0,-i 

S) G 23 ( 


_ i I HC_ | 

9 ' K 


2 « r (| + ^ 


+ £f 1 


0 


1-4 


V* 




ijL) G 23 ( 0 + -|+£f I 1 ) 


(63) 


( 


1 - 2 


2 « r (| + 3f) 


G 23 fe-i+^l 1 )) 


Therefore we see that the vertical component (y t ) of the tip of the SLE trace shows serious deviations 
from the ordinary SLE results, in contrast to the horizontal component which has the same dependence 
(with modified proportionality constant): 


y t ~ v / 2 Kt + ^V^Kt 


x t = fy/2Kt 

-1 


(64) 


4.3 CFT Prediction for the PDF of SLE(tt, p c ) Traces 


As SEC 3.2 in this section we apply the CFT constraints to the correlation function which are repre¬ 


sentative of the probability density of the SLE(k, p c ) random curves. For the this case we can use the 
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Eqs [35] and [34] to obtain 


[j3c„ “ 2/ ^ Re 


1 


1 d 


1 d 


1 


a; 00 -Co 


(z - Co) 2 J z - Co dz z- Co <9z 

]P(a;,y,Co,a; 00 ) = 0 . 




(65) 


After some straightforward calculations, using the fact that P = 2 \ F h{u) (Eq. 36) we obtain: 


k 2 , 2 (z - Co) fl 

2 Co (x- Co) 2 + y 2 


2 V 




6 — K 


( X - Co) 2 + 2/ 2 " v ' X°° - Co " c ° ' 2:00 - Co 


d :r 


8 h 


F {x 2 + y 2 ) 2 


h(u ) = 0 
(66) 

Rearranging the derivatives in terms of u = ^ = ^ + x°°-Co ’ aRer some algebra 

it is concluded that: 


^ (l + it 2 ) 2 <9 2 + 4u (l + u 2 ) d u + 8/if /i(it) = 0 


( 67 ) 


which is the same as Eq w8] by replacing y —> u. Therefore for the case SLE(k, p c = k — 6) we have 


P(x,y,Co,x°°,e) = (jj) M u ) = (f) f (l + u 2 ) ( ® h)/(2h) . 

This result is expected, since due to conformal invariance of the model we can generate SLE( k, k — 6) 


curves from chordal ones by the map z = <j){w ) = X cx,™ w in which z indicates the coordinate in which 
the chordal SLE grows and w = x + iy the dipolar one. Therefore we expect that ( e’ is the transformed 
radius under <j>) 


Pdipolar('a , 5 X ) — -/"chorda I (0(uj) ) — 


We can easily verify that = — u and since h(u) is even with respect to u we find that 

/ //\ 2 df f x°° )2 -|-y2 A 

Pdi P oiar(»i x°°) = I y J h(u) in which e" = -- ^2 --e 7 is the renormalized radius. 


Vlm[0(ty)] 


f h(r J = 

[1 Im[0(«7)] ’ 


( 68 ) 


5 Conclusion 

In this paper, we analyzed the statistics of the SLE curves described by SLE(/t, p c ). We obtained the 
Fokker-Planck equation governing these SLE traces. We found that the long-time limit behavior of 
the x and y components of the tip of the growing curves are not the same. The x component is only 
renormalized (meaning that the proportionality constant changes), and the y component is changed 
considerably, i.e. its functional form changes. The exact forms of these changes were reported in this 
paper. 
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